V. V. KOZLOV 

Steklov Institute of Mathematics 
Russian Academy of Sciences 
117966 Moscow, Russia 
E-mail: kozlov@pran.ru 



7?&C 



NOTES ON DIFFUSION IN COLLISIONLESS MEDIUM 

O . Received October 06, 2003 DOI: 10.1070/RD2004v009n01ABEH000262 



On 



> 



> 



A collisionless continuous medium in Euclidean space is discussed, i. e. a continuum of free particles moving inertially, 
without interacting with each other. It is shown that the distribution density of such medium is weakly converging to 
zero as time increases indefinitely. In the case of Maxwell's velocity distribution of particles, this density satisfies the 
well-known diffusion equation, the diffusion coefficient increasing linearly with time. 
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1. Dynamics of collisionless medium in the Euclidean space 



We are going to consider a very simple object — a collisionless continuous medium, i. e. a continuum 
of free particles moving inertially, without interacting with each other. The configurational space of 
a particle is the Ti-dimensional Euclidean space R n with orthogonal coordinates xi, . . . , x n ; let R n = 
, = {o>i, . . . ,u! n } be the velocity space. The direct product R" x R™ = T is the phase space of a free 
particle. 

' Let p(u>, x) be the particle distribution density at the initial time t = 0. We can assume that the 

density is normalized to the total mass of the collection of particles (the total mass is supposed to be 



in 
o 



e ■ f 



finite). In other words, p is the density of some probability measure: p ^ and 

pd n ud n x = 1. 



From the very beginning, it is possible to use the probabilistic approach and treat a collisionless 



space i« n 



^ ■ medium as a Gibbs ensemble of identical systems, where each system is a free particle in the Euclidean 
je R n . 

According to the elementary principles of statistical mechanics, density pt at time t is given by 

Pt(oj, x) = p(u>, x — Lot). 
It is clear that po = p. It is also clear that 

u(x,t) = / p(oj, x — ujt)d n uj (1-1) 



is the density of the collisionless medium at point x at time t. Our objective is to study the concen- 
tration u, its evolution and limit behavior as t — ► dboo. 

In Ref. [1] , the problem of evolution of a collisionless medium inside a box with mirror walls was 
discussed. Upon a simple regularization, this problem can be reduced to the problem with periodic 
boundary conditions: the configurational space R n = {x} is factorized using the lattice (2-7rZ) n . As a 
result, the phase space T is the direct product of torus T n = {x±, . . . , x n mod2ir} and R = {u}. 
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Theorem 1. Let p G L\{T) and (p be the characteristic function of a bounded measurable re- 
gion D C R n = {x}. Then 

J u(x,t)(p(x)d n x = J u(x, t)d n x (1.2) 

B™ D 

as t — > ±oo. 

This theorem is intuitively obvious: the particles scatter to infinity, each with its own velocity, 
and, therefore, their concentration in any finite region of K n = {x} is decreasing indefinitely. 

In fact, (1.2) holds for any essentially bounded measurable function <p : lR n — > R. Equation (1.2) 
has the following meaning: weak density limit pt is zero as t — ► ±oo. This assertion was proven in 
Ref. [1] for a compact configurational space. In any case, 

lim / p(u, x - Lut)ip(x)d n xd n uj = / p~(pd n xd n uj, (1.3) 

t^±oo J J r 

r 

where p is the Birkhoff average of p: 

T 

p(u,x) = lim / p(u,x — ut)dt. 
t^oo It J 

— T 

It is easy to calculate: 

X+U)T 



p=\im 1 ±- [ p(w,Ode = 

T— >oo ZWT 7 



for almost all u, because (according to Fubini's theorem) the integral of p(uj, £) over the variable 
exists (and is finite) for almost all uj. Specifically, the integral (1.3) is also zero. 

The question whether u(x,t) itself tends to zero as t — > ±oo is somewhat more interesting. We 
discuss it for the case where p is the product of two summable functions h(uS) and (f(x). Then, the 
integral (1.1) takes the form 

u(x,t) = / h(u)tp(x - cut)d n uj. (1.4) 



Let us present <p in terms of Fourier transformation 



and put 

H(z) = _L I h{u)e- l ^ uj) d n oj. 



(V27r) n Jr™ 

Theorem 2. Suppose that & is a bounded summable function and H G Li(R n ). Then, for 
all x G R n , function (1.4) tends to zero as t — > dboo. 

Indeed, according to Fubini's theorem, 

L J $(£) e *(^) J h(io)e- it{ul ^d n LO d n £ = J $(£).H'(tf)e i ( x '* ) d n f. 



When t > 0, this integral is equal to 



i J H(z)${^)e i(x > z)/t d n z, 
i" 

which is of order 0(1 /t) as t — ► oo (by the assumption of the theorem). 
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2. Heat conduction equation 

The equation is 



u T = a 2 Au, (2.1) 



where r is a time variable, A is the Laplace operator, while a = const. Equation (2.1) is a special 
form of the diffusion equation; a 2 is the diffusion coefficient. 
The solution of (2.1) is well-known 

\*-Z\ 2 

u(x, t) = (2tj J_ )n J e ^ <p(Z)d n Z, r > 0, (2.2) 

r 

where up is the initial temperature distribution, and |g| = q\ + . . . + q 2 . The function p is customarily 
supposed to be continuous and bounded. The latter condition ensures convergence of integral (2.2), 
while the continuity property allows proving that 

lim u(x, z) = ip(x). 

The point is that the exponential term in (2.2) (together with the term outside the integral) tends to 
the delta- function 5(x — £) as r — > 0. 

However, the integral (2.2) also converges on the assumption of summability of op (i. e. when 
ip G L\). It turns out that (2.2) can be presented in the form of (1.4), and this, among other things, 
implies that u(x, 0) = p(x) if p is summable. 

Let r = t 2 /2 and ^ = x — lot, where to = {uj\, . . . , u> n ). Then, d£j = —tdu>j and 



|g-ei 

/ - 2a 2 t 



(-iy 



p(£)d n t 



-oo |„,|2 



(aV27r) n 



oc 
|u;' 2 

1 



(V27r<r) r 



e 2 <T 2 <p( x - ujt)d n uj = (2.3) 

2 ^ 2 <p{ X - UJt)d n LV. 



In particular, u(x, 0) = <y?(x). 

Thus, if we adopt the normal law of velocity distribution 



.12 



1 



/i(w) = — = e 2 - 2 , (2.4) 



(v 27rcj) 

then the density u(x,t) of the collisionless medium, given by the integral (1.4), satisfies the diffusion 
equation 

ut = ta 2 Au. (2.5) 

The diffusion coefficient to 2 increases indefinitely with time. As distinct from the heat conduction 
equation, this equation is invariant under time reversion t \-> —t. This reflects the property of re- 
versibility of the equations of motion for a free particle. Specifically, concentration of particles at any 
point x G M. n decreases indefinitely both as t — > +oo and as t — > — oo. 

From the statistical mechanics point of view, it would be more appropriate to treat the distri- 
bution (2.4) as a Maxwell distribution, dispersion a 2 being proportional to the absolute temperature 
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of the gas. This distribution does not vary with time (because the medium is collisionless), and the 
temperature field is proportional to the density of the collisionless medium (after identifying t 2 /2 
with r). 

The simple equation (2.3) is also useful for the analysis of heat propagation in R n . For example, 
let (p > inside an open bounded region D C W 1 and tp = outside this region. Then u(x, t) > at 
any point for arbitrarily small t > 0. This property immediately follows from the law of motion 

of a collisionless medium : however distant a point x G R™ may be, it will be reached in an arbitrarily 
small time by very fast particles, located initially in D. 

3. An example of a nonstandard diffusion equation 

It would be a mistake to think that functions of the form (1.4) satisfy the diffusion equation in its 
well-known form. Let us put, for example, n = 1 and 

h(uj) = when to ^ and h(to) = when uj > 0. 

Then (1.4) becomes 

o 

u(x,t) = J e w <p(x — ust) dw . (3-1) 

— oo 

Integrating by parts yields the equation 

u = tu x + (p(x), (3-2) 

which does not contain the derivative of Ut at all. Putting t = in (3.2), we find that u(x,0) = <p(x). 
However, Theorem 2 cannot be applied straightforwardly because 



H(z) 



1 



y/2n J v / 2tt(1 - iz) 



does not belong to Li(IR). However, if, instead of the boundedness of |$(z)|, we require that this 
function vanishes at infinity as 0(|z| _Q! ), a > 0, then again we can say that the integral (3.1) tends to 
zero as t — > ±oo, for each i£K, 

The function tp can be eliminated from (3.2) if we replace (3.2) by the equation 

ut = (tu x ) t (3.3) 

and the Cauchy condition u\t=o = <p(x). Thus, (3.3) should also be considered a diffusion equation. 

Equation (3.3) is not invariant under time reversal. This fact can be easily explained: all the 
particles of a collisionless medium move to the left. To have symmetry between the past and the 
future, one should assume that there is symmetry between "left" and "right" in the distribution of 
velocities. We come to a simple 

Proposition. If h(—uS) = h(u>), then u(x,—t) = u(x,t). 

Indeed. 



u(x,—t)= I h(uj)ip(x + ujt) dui = — J h{—u)')ip{x — u't) dw' 

oo 

h{uj)ip{x — tot) dt = u(x, t). 



-oo 

oo 
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h{u>) = \e-\ u \. 



Integrating by parts, we obtain 
o 

e"V(x - tot) duo = <p(x) + t<p'(x) + t 2 I e"up"{x - uot) duo, 



— oo — oo 

oo oo 



e w up(x — tot) duo = tp(x) — tup'ix) + t 2 J e ^up" {x — uot) duo. 





This leads to 



t 2 

u = up(x) + -ju xx , 



t 



with u\t=o = f(x). This equation is equivalent to the equation of evolution 

2 

with the initial condition u\t=o = (p(x). Unlike (3.3), this equation is invariant under time reversal. 

4. Diffusion in the compact case 

Now, let the configurational space be an n-dimensional torus T™ = {x\, . . . , x n mod 2tt}. In this case, 
the density of a collisionless medium is also given by (1.1). As it was shown in [1], the weak limit 
of u(x, t) as t — > ±oo is equal to 



p(uo,x)d n uod n x. (4.1) 



Under certain additional conditions, we can state that u(x,t) —> u as t — > ±oo for all x E R n . To 
this end, suppose that the density u(x,t) is given by the integral (1.4). Let 

J2^me i{m ' x) , meZ" (4.2) 

be the Fourier series of a bounded measurable function up: T n — > R. 
Theorem 3. Suppose that the series 

converges; then u(x, t) — > u as t — > ±oo for all x G T n . 
Inserting (4.2) into (1.4), we obtain: 

u(x, t) = u+J2 W i{m > x) J /i(w)e- i(m '^ d n uo, (4.4) 

m^O R n 

where 

u = (fo h(uo) dTuo 
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coincides with (4.1). Because of the convergence of the series (4.3) the summation over m and inte- 
gration over u can be interchanged. 

Let e be an arbitrary small positive number. Since the series (4.3) converges, there exists a 
number N, depending on e, such that 



Vme i{m,x) J /iHe _i(m,£j) * d n uo 



|m|>iV 




h(u) d n uj 



is less than e/2. Then, the terms with indices subject to \m\ ^ N and m/0 tend to zero as t — > oo, 
because (by the Riemann-Lebesgue theorem) so does the integral 

/i(u;)e- i(m ' w) * d n to 



Thus, there exists T{e) such that for t > T the sum of this finite number of terms is less than e/2, 
which is the required result. 
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